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XXXII. On the Differential Coefficients and Determinants of Lines, and their Application 
to Analytical Mechanics. By A. Cohen, Esq. Communicated iy Professor Stokes, 
Sec. U.S. 

Eeeeived May 8,— Bead June 19, 1862. 

Chapter I. 

1. I propose in these pages to prove the principal theorems of dynamics in a manner 
which appears to me both simpler and more methodical than that in which they are 
generally proved; and I believe that I shall be able, by applying a few conceptions 
which spring naturally from the principles of higher algebra and statics,' to give a clear 
interpretation to most of the more complicated formulae in dynamics, as well as to the 
several analytical steps which lead to those formulae. 

2. There are many reasons why the diagonal AD, which is constructed on the straight 
lines AB, AC, should be considered as the sum of those two lines. Those reasons may 
be found developed in De Morgan's ' Double Algebra,' in Warren ' On Imaginary 
Quantities,' and in the Tract of Benjamin Gompertz ' On Imaginary Quantities.' . 

I shall therefore call AD (AD being the diagonal of the parallelogram constructed on 
AB and AC) the complete sum of AB and AC, and the two lines AB and AC will be 
called the components of AD. Moreover, denoting AB, AC, AD by P, Q, K respectively, 
I shall express their relation to one another by the equation 

K=(P)+(Q). 

I shall likewise denote by (— Q) a line equal and opposite to Q, and define (P)— (Q) 
to be the same as (P)+(—Q), calling it the complete difference of P and Q. 

All lines which have the same length and direction will be considered as equal to one 
another, so that any line is equivalent to a line through the origin having the same 
length and direction. 

3. It evidently follows from the above definitions, that the complete sum of AB and 
BC is AC, and the complete difference of AB and AC is BC. 

4. Suppose now Q to represent a line which varies with the time t both in length and 
direction. The complete difference of the two consecutive values of Q after an increment 
of time At may be called the complete increment of Q, and may be denoted by A(Q). 
Moreover, if we divide the length of A(Q) by At, and take the limit of that ratio, then 
the line which has that limit for its length, and which has for its direction the direction 
of A(Q), when At diminishes without limit, will be called the complete differential 
coefficient of Q, and will be denoted by D/Q). 

5. It will be sometimes found convenient to denote a line of length r, which is parallel 
or perpendicular to a line P, by (r || to P) or (r_[_ r to P). Moreover, if n represent a 
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470 ME. A. COHEN ON THE DIFFEEENTIAL COEFFICIENTS 

numerical quantity, then nP may be used to denote a line which is in the direction of 
P, and whose length bears to that of P the ratio of n to 1. 

6. The following Lemmas, which will be of constant use, are all but self-evident : — 

I. If K=(P)-(Q), then (R)-(P)+(Q)=0. In short, the ordinary rule of signs holds 
good. 

II. (nT)±(nQ)-n{(P)±(Q)\. 

III. The projection on any line or plane of the complete sum or difference of two 
lines is equal to the sum or difference of their respective projections on the line or plane. 

7. Whenever there is no risk of any mistake, the brackets may be omitted in the above 
and similar formulae, and the complete sum or difference of lines may be spoken of simply 
as their sum or difference. 

What has been hitherto said may be of course extended to all magnitudes whatsoever 
which can be adequately represented by straight lines, such as forces, velocities, axes of 
couples, axes of rotation, and accelerations, &c. 

8. The application which can be made in dynamics of this conception of the 
complete differential coefficient of a line, will become at once apparent from the follow- 
ing considerations. 

Suppose a particle to be moving from A to B. Let O be any fixed point. Then the 
particle's velocity is, according to its very definition, represented in 

AB 

magnitude and direction by the limit of -rr. But AB is the com- 
plete difference of O B and O A, or the complete increment of the 

AT? 

radius vector O A, and therefore the velocity, being the limit of -gj, is the complete 

differential coefficient of the radius vector. 

In the next place let A and O B in the last figure represent in magnitude and direc- 
tion the successive velocities of a particle at times t and t-\- At respectively. Then, since 
OB=(OA)+(AB), it follows that, if with the velocity OA we compound the velocity 
AB, we shall obtain the velocity OB, and therefore the particle's acceleration is repre- 

AB 

sented by the limit of -tt, or by the complete differential coefficient of the velocity OA. 

Hence we have the following proposition : — 

If a particle's velocity and acceleration be represented by straight lines, the velocity 
will be represented by the complete differential coefficient of the radius vector drawn 
from a fixed point to the particle, and the acceleration will be represented by the complete 
differential coefficient of the velocity. Or more briefly, the velocity is the complete dif- 
ferential coefficient of the radius vector, and the acceleration is the complete differential 
coefficient of the velocity, and is therefore the second complete differential coefficient of the 
radius vector. So that if E, denote the radius vector, and V and F denote respectively 
the velocity and acceleration, we have 

V=tD,(E), F=D,(V)=D*(R). 
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9. Such then being the connexion which exists between the differential coefficient ot 
the radius vector and the velocity and acceleration of a particle, I will proceed to prove 
some of the principal propositions concerning the differential coefficients of lines, and to 
apply them to the dynamics, or rather the kinematics of a particle. 

The first proposition is the following : — 

If P, Q, E represent straight lines, and if we have 

(P)±(Q)=E, 
then 

D,(P)+D,(Q)=D/E). 

For let P, Q, B after an interval of time At become P', Q', B/ respectively, then we 

have 

(F)+(Q')=E', 

and therefore, by Lemma I. of section 6, it follows that 

P')-(P)} + KQ')-(Q)} =(£')-(*), 

or 

A(P)+A(Q)=A(E). 

Therefore by Lemma II. of section 6, we have 

/A(P)\ /A(Q)\_A(R), 
\ A* / ^ ^ A* / — A* ' 

and taking the limit of both sides of this equation, we obtain 

D t (P)±D,(Q)=L\(E). 
Similarly it may be shown that 

D < {(P)±(Q)±(R)}=D,(P)+D,(Q)±D ( (E). 
Moreover, denoting the second complete differential coefficient by D,, it follows that 

D?{(P)±(Q)±(R)}=D i {D,(P)±D,(Q)+D < (E)} 
=D?(P)+D?(Q)±D?(E). 

10. Suppose now a line Q to have Q,,,, Q^, Q,. for its components parallel to the axes 
of coordinates Ox, Oy, Oz ; it is evident from Lemma III. of section 6, that 

Q=(Q.)+(Qr)+(a). 

It follows, therefore, from the preceding section, that 

D,(Q)=D < (Q,)+D < (Q,)+D t (Q,) (I.) 

and 

D?(Q)=D|(Q.)+D?(QJ+D?(Q,). (II.) 

These equations are true whether the axes of coordinates are fixed or move. But 
supposing the axes to he fixed axes, let q m q p q z he the respective lengths of Q„ Q^,, Q*. 
Then it is evident that, as the direction of Q, does not vary, D^Q,) is a line whose 

direction is that of Q* or Ox, and whose length is -j£ ; and similarly, D^Q*) is a line whose 

at 

direction is that of Ox, and whose magnitude is -~. Similar results hold good for 
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472 ME. A. COHEN ON THE DIFFEKENTIAL COEFFICIENTS 

Di(Qy), D?(Q^,), &c. Therefore the two equations (I.) and (II.) evidently show that the 
components of T> t (Q) and D*(Q) parallel to Ox are respectively equal to -J? and -~- 

11. It is easy to apply the above results to the velocity and acceleration of a particle. 
For let Q in the last section stand for the radius vector of a moving particle, then the 
components of the radius vector are respectively equal to x, y, and z ; and since the 
velocity is the complete differential coefficient, and the acceleration is the second complete 
differential coefficient of the radius vector, it follows from the last section that the com- 
ponents parallel to Ox of the velocity and of the acceleration are respectively equal to 

ttl* fl if* /a 

-in and -is-. So that if v M v y , v x be the components of the velocity, and f m f y , f K be the 
components of the acceleration, we have the elementary formulae 

St OfU (iz 

%=j- p and similarly «,= #, ««=# \ 

n d*x t • m n /> cPy „ d^z 
f*=jp> and fiamilarly/,=^,/.=^ 

12. Our next proposition will arise from investigating the complete differential coeffi- 
cient of a line Q, which varies both in magnitude and direction with the time t. 

Let Q at time t be the line O A, and let it become O B at time 

t-{-At, so that we have 

A(Q)=(OB)-(OA)=AB. 

Produce OA to C, making OC=OB, and draw BC. Let OA 

and OB have for their respective lengths q and q-\-A([, and let angle B O A =05. 

Then 

A(Q)=AB=(AC)+(CB). 

Therefore, by Lemma II. of section 6, 

A(Q)_/AC\ /CB\ 

-&-=\At) + (jt)- ■■•■■■■■■■ (I.) 

Now diminish At indefinitely and take the limit of the last equation. The limit of 

A(Q) AC 

. , is L\(Q), the complete differential coefficient of Q. The limit of -rr is evidently 

AC da 
a line whose length is the limit of -tj, or -£, and whose direction is that of OA or of Q. 

CB 

Finally, since COB is an isosceles triangle, -rj has for its limit a line whose length is 

OA limit of -7-:, and whose direction is perpendicular to OA or Q, and in the plane in 
which Q is moving at time t ; and if rs be the rate at which the direction of Q is varying 
at time t, w=limit of 7-7. Therefore, taking the limit of equation (I.), we obtain 

D,(Q)= (§ || to Q) +(£m±* to Q), 
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the latter line (qvs _|_ r to Q) being in the plane in which Q is moving at time t. This is 
the fundamental proposition concerning the differential coefficient of a line, and may be 
stated in the following form : — 

The complete differential coefficient of a line Q, whose length is Q and whose direction 
is at time t varying with an angular velocity vx, is the complete sum or is compounded of 

two lines, one being -n in the direction of Q, and the other being qyr in a direction perpen- 
dicular to Q and in the plane in which Q is moving at time t. The former of these two 
lines would evidently be the complete differential coefficient of Q, if the length of Q only 
varied, and the latter would be its complete differential coemcient if the direction of Q 
only varied ; and in this sense it may therefore be said that the complete differential 
coefficient of a line is the complete sum of the two partial differential coefficients obtained 
by varying separately the length and the direction of Q. One of these partial differential 
coefficients may be called the length-differential coefficient, and the other the direction- 
differential coefficient of Q, and the complete sum of these two constitutes the complete 
differential coefficient of Q. 

13. Let Q in the preceding section stand for the velocity of a moving particle. Then 
D t {Q) will be the particle's acceleration, q will be the velocity v, and the direction of Q 
will be that of the tangent to the particle's path. Finally, rsdt will be the angle between 

ds 
two consecutive tangents, so that nrdt=—, ds being an element of the particle's path, 

1 uS V 

and g> the absolute radius of curvature. Therefore &=--}:=-. It follows then at once 

dv 
from the last section, that D,(Q), the particle's acceleration, is compounded of -r, along the 

tangent, and vzs- or — perpendicular to the tangent and in the plane in which the radius 

vector is moving at time t. In other words, the resolved part of the acceleration along 

dv d s v 

the tangent is Tf = Jfi , > an( l the resolved part along the absolute radius of curvature is — . 

14. The same fundamental proposition of section 12 enables us to investigate D,(Q), 
the second complete differential coefficient of a line Q, if we suppose that line to move 
always in the same plane. We have, namely, 

D,(Q)= (J 1| to Q) +(ms±' t0 Q)- 
Now, in order to find D,(Q), we must ascertain the complete differential coefficients of 
( -J || to Q ) and of (™aj_* to Q). The complete differential coefficient of the former line 

(-£ || to Q) is, by section 12, 

(Sl|toQ)+«(JrtoQ). 
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Again, the complete differential coefficient of the line {qys JJ to Q) is similarly the 

complete sum of the line (-j t (^)jJto QJ and of a line whose length is q&* and whose 

direction is perpendicular to the line (§w X r to Q), and whose direction is therefore 
evidently opposite to that of Q. Hence Df(Q) is the complete sum of 

^f |[ to Qj and (w-f X r t° Q) and (j t (q™)j_'to Q) and (— q^ || to Q). 



Therefore 



D ?(Q)={(5-^ 2 ) II toQ}+{(s*+3 w) r toQJ- 



In other words, the components of DJ(Q) parallel and perpendicular to Q are respect- 
ively 

15. This last result may be easily applied to the dynamics of a particle. For, let Q 
stand for the radius vector of a particle moving in a given plane. Let that radius vector 
have r for its length and w for its angular velocity ; then, since the acceleration equals 
the second differential coefficient of the radius vector, it follows at once from the last 
section, that the components of the acceleration parallel and perpendicular to the radius 
vector are respectively 

l?- rrs and r7g +Jt H» or - Jf (?■%). 

16. This last result is, however, but a particular instance of the connexion which 
exists between the actual motion of a particle, and its motion relatively to axes which 
move in the same plane as the particle moves. It will be found that that connexion 
may be easily deduced from the solution of the following problem : — 

" Supposing the axes of coordinates Ox and Oy to move about O in the plane of xy 
with an angular velocity w at time t, it is required to find the complete differential coeffi- 
cient of a line Q which moves in that plane, the lengths of Q's components along the 
moving axes being given." 

Let Q have for its components Q, and Q^,, and let the respective lengths of these be 
q„ and q r Then, by Lemma III. of section 6, we have 

Q=(CU+(QJ. 
Whence it follows that 

D,(Q)=L\(Q,)+D,(Q,). 

Now since Q, and Q^ vary in direction as well as magnitude, and since the angular 
velocity of their change of direction is &, we have, by the fundamental proposition in 
section 12, 

D '( Q *)= (% II t0 0a ) +( w & -i-* t0 °*)» 



D ^=(ll H to %)+( w ^J- r to Oy). 



(I-) 
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But since the lines (tsq t j_ r to Ox) and (ra^ j_ r to Oy) are respectively proportional and 
perpendicular to Q^ and Q^, and since the latter lines have Q for their complete sum, it 
evidently follows that the former two lines have for their complete sum a line which is 
perpendicular to Q, and whose length is &q. 

Hence 

r> ^ = (^ll tG0 ^) + (^ll to0 2')+^ J L r toQ).. . . . (II.) 

17. This last formula is true whether the axes be rectangular or oblique, and may be 
made the basis of all the formulae of relative motion in one plane. 
It may be observed that the line 



(tll'° *) + Ct» to °2') 



is what would be the complete differential coefficient of Q if the coordinate axes were 
fixed ; and it may therefore be called the complete differential coefficient relative to the 
moving axes, or, more briefly, the relative differential coefficient of Q. So that the above 
formula shows that the complete differential coefficient of Q is its relative differential 
coefficient together with a line {■&% J_ r to Q), the latter line being drawn towards the direc- 
tion in which the axes are revolving. 

18. If the axes of coordinates be rectangular, then the line (rsq,, J_ r to Ox) is evidently 
the same as (mq x || to Oy), and the line (bj<^ _)J to Oy) is the same as (— mq 9 1| to Ox) ; and 
therefore, looking at the equations (I.) in section 16, we see that 
D,(Q)=D,(Q,)+D,(Q,) 



:((!-.?,) || toO*) + ((!»+. fc ) || toO, 



In other words, the components of L\(Q) parallel to 0* and Oy are respectively 

I— * - d %+n- • • Pit) 

The same result may be also deduced from formula (II.) in the same section, if we 
resolve the line (mq JJ to Q) along the rectangular axes of x and y. 

19. Let us apply the above formulas first to the velocity of a particle. 

Suppose, then, a particle to move in a given plane, and that the rectangular axes of 
coordinates in that plane revolve about the origin with an angular velocity & at time t. 
Let v„ and v y be the components of the particle's velocity along the moving axes. Then, 
since the velocity is the complete differential coefficient of the radius vector, and since 
x and y are the components of that radius vector, it follows from the formulae (III.) of 
the last section, that 

dx 

dy , 
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If the radius vector be chosen as axis of x, then x=r, y=0 ; therefore % == ~it along 

the radius vector, v y =r i m perpendicular to the radius vector, where nr is the angular 
velocity of the radius vector. 

20. Let us now apply the same formulae (III.) of section 18 to the acceleration of a 
particle. Let v x and v s , as before, denote the components of the velocity, and let f x and f g 
denote the components of the acceleration of the particle. Then, since the acceleration 
is the complete differential coefficient of the velocity which has v„ v y for its components, 
it follows at once from the formulae (III.), that 

j} do 

Suppose now the axis of x to be the radius vector, then we have already shown that 
v x =-n, v y =rw. Therefore by substituting these values in the last formulae, we see that 

d?r 
f„ the acceleration along the radius vector, is -p — rvs 2 , and f y , the acceleration perpen- 
dicular to the radius vector, is -i (m)-\-is -t=--t (tot 2 ), which is the same result as was 

obtained in section 15. 

21. Eeturning to the more general case, we have, as before, 



dt, 

dVy 
It' 






and substituting in these the values already obtained for % and v y , namely, -jj.—ym, 



-tt-^-Xis, we find 



j, (fix n dy . d'as 

j, d*y dx dvr 

fv=-dw+^-dt™-™y+ x Tv 



(IV.) 



which are the formulae for the components of the acceleration in terms of the coordinates 
of the particle. 

22. The last formulae (IV.) may also be obtained in the following manner. 

Let Q^, and Q y represent respectively the components of the radius vector Q along the 
axes of x and y. Then 

Q=(Q.)+(Q,). 

Therefore 

D 2 t (Q)=D?(cg+D 2 (q,). 
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But since the axes of x and y revolve at time t with an angular velocity w, and since 
Q # and Q r have x and y for their respective lengths, it follows from section 14, that 



and that 



D?(Q,)={(§-^ 2 ) || toO*}+[^ (aV)rtoOr}, 



Whence it is easy to see that the components of D?(Q), or of the particle's acceleration, 
are 



j, d*x 1 d . „, 



m 



which equations are clearly the same as those obtained in the preceding section. 

We shall soon prove similar formulae for the more general case of a particle and axes 
of coordinates moving in any manner whatsoever in space of three dimensions, and 
therefore, in order to prevent needless repetition, we shall postpone the further discussion 
and complete interpretation of the equations (IV.) or (V.). 

23. It is, however, interesting here to observe that all the results already obtained may 
be readily deduced from the principles of what Professor De Morgan has called "Double 
Algebra." According to those principles, namely, the radius vector B whose length is 
r and whose inclination to a fixed line is 6, is symbolically represented by ri ^ _1 , so that 
we have 

Therefore 

and the last expression represents the complete sum of 



(J || to K) and (rgrtoB). 



This result is the same as that arrived at in section 12 
Again, D?(B)=D,(D,(B))= 

v-r /d*r 



t + 



-/d s r , ,— -=- d / dt\ , dt ,—= fdr . dh /— tAN 
1 {dfi+^^di^diJ+Jt^- 1 {Jt+ r dt^- 1 )) 

.,,— /<Pr df , ,— - fd f dt\ dra\\ 



and this expression represents the complete sum of a line 

<2V dt n , ~ . ,. d f di\ . drdd 



dt*~ r Tt Ii toR and a lme di{ r dt)+Jtdt-i- ton - 



This result is the same as that arrived at in section 14. 
mdccclxii. 3 T 
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Finally, in order to obtain the formulae for relative motion, we have merely to put 

where & is the angle made by R with the moving axis of x, and a is the angle made by 
that moving axis with a fixed line. It follows then that 

D,(R)=s» '-'!(« ^)+ r ^y/-l ,<»+«W=i. 
Now it is evident that s* ^~ 1 j \n e ^~ 1 ) represents the relative differential coefficient of 

R, and r^-\/ ~*l s(9+ " represents a line r^ J_ r to R. We thus obtain the same 

result as in section 16. 

By differentiating again it would also be easy to deduce the result of section 20, if we 

observe that i a ' J ~ 1 ^\rs e ^~ 1 ) represents the particle's relative acceleration whose corn- 
ea? , d% 
ponents are -™ and -rf • 

Chaptee II. 

24. In order to extend the formulae which we have proved for the motion of a particle 
in one plane to the motion of a particle in space, it will be found very convenient to 
make use of a conception which presents itself in statics, as soon as the equilibrium of 
a solid body is treated of in that science. 

Let O A and O B be any two straight lines drawn from the origin O. ^^ 
If then O A represent a force P, and if we apply at B a force — P, we 
shall obtain a cowple. Let O C be the axis of that couple. We know 
then from statics that, if O A and O B have for their projections on the 
axes of coordinates X, Y, Z and x, y, z, then O C has for its projections 

zY-yZ, xZ—zX, yX-xY. (I.) 

Now the relation which the line O C bears to the lines O A and O B is one which not 
only presents itself in statics, but which also plays a very important part in the differ- 
entiation of lines, and in the dynamics both of a particle and of a body. For this reason 
it will be proper to treat of the relation in question quite independently of statical 
considerations ; and since the expressions (I.), which are the projections of O C, are evi- 
dently what are called determinants, I shall call the line O C the determinant of B to 
OA. 

Hence we have the following definition : — 

" The determinant of a line Q to a line P is a line which is equal to twice the area of 
the triangle of which the lines P and Q drawn from the origin are sides, and which is 
perpendicular to that area, and the line is moreover drawn in such a direction that, to 



-p 
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an eye looking along it towards the origin, the revolution of Q towards P appears to be 
a revolution in the positive direction." 

The determinant of Q to P may be briefly denoted by 

det (Q, P). 

It is evident from the above definition, that the determinant of Q to P is a line equal 
and opposite to the determinant of P to Q. 

Moreover, if the projections of P on the axes of coordinates he p M p v , p„ and those of 
Q be q m q v , q z , then it follows from the formulae (I.), that the determinant of Q to P or 
det (Q, P) has for its projections 

g&.—2ffi r ft?.— &P« 2&—2*P* ( IL ) 

25. The connexion which exists between the notion of a determinant of lines and the 
elementary conceptions of dynamics may be easily made apparent. For suppose a 
particle at the extremity B of O B to be revolving about the line O A with an angular 
velocity represented in magnitude by O A, then if O C be drawn perpendicular to the 
plane A O B, and equal to twice the area A O B, it is evident that O C will represent 
the linear velocity of the particle. But O C is then by definition the same thing as the 
determinant of O A to O B. Whence it follows that the determinant of O A to O B 
represents the velocity of the point B, due to a rotation whose axis and angular velocity 
are represented by O A. 

This result, together with the result of the preceding section, may then be recapitulated 
in the following manner. If V represent the velocity of a particle at the extremity of 
the radius vector R, and the particle rotates about an axis which is represented by the 
line n, then, if the angular velocity is represented by the length of ft, we have 

V=det (H, R). 

Secondly, if P represent a force at the origin and R represent the radius vector at the 
extremity of which a force — P acts, then the axis of the couple (P, — P) is det (P, R) 
or det (R, — P) ; so that det (R, P) is what French writers call " the moment-axis of a 
force P with respect to the origin." 

26. Such, then, being the connexion between determinants and statical and dynamical 
conceptions, I will proceed to prove some of the more important propositions concerning 
the determinants of lines. 

The most important theorem concerning the determinants of lines is the following : — 
" If P, P', and Q be three straight lines drawn from the origin, then 

det (P, Q)+det (P, Q)=det {(P)+(P), Q}." 
This proposition might be easily proved by geometry, but it is at once deducible from 
statics. For, consider two couples having a common arm Q, and having forces P and P' 
respectively acting at the extremity of Q at the origin O. The resultant of those two 
couples will be a couple having the same arm Q, and having for its force acting at O the 
resultant of P and P', or (P)+(P'). Now it is proved in statics that the axis of this 

3t2 
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resultant couple is the complete sum of the axes of two component couples. Therefore, 
substituting for those axes the equivalent determinants, we see that the determinant of 
(P)-f-(P') to Q is the complete sum of the determinant of P to Q, and of the deter- 
minant of P' to Q. Thus we have 

det(P,Q)+det(F,Q)=det{(P)+(F), Q} (I.) 

And it may similarly be shown that 

det(P, Q)-det(P'Q)=det{(P)-(F), Q}. . ........ (II.) 

The same proposition follows also easily from a consideration of the linear expressions 
(II.) in section 24 for the projections of a determinant, and is in fact equivalent to a 
fundamental theorem concerning algebraic determinants, which theorem may be found 
in Salmon's ' Lessons on Higher Algebra,' section 19, page 9. 

27. The proposition proved in the last section will be found to be of constant use in 
explaining and shortening analytical processes in mechanics. One useful application can 
be made of it in proving " the parallelogram of angular velocities." For taking Q to be 
the radius vector of a particle, and P and P' to represent two axes and angular velocities 
of rotation then the formula (I.) of the last section translates itself at once by means of 
section 25 into the following proposition : — " The linear velocity of a particle due to a 
rotation whose axis and angular velocity are represented by the line P, compounded 
with the linear velocity due to a rotation similarly represented by the line P', is equiva- 
lent to the linear velocity due to a rotation represented by the complete sum of P and 
P'." And this is evidently the same as the proposition called " the parallelogram of 
angular velocities of rotation." 

28. Let us next investigate the complete differential coefficient of det (P, Q). 

We will first premise that, if m be any numerical quantity, it follows evidently from 
the definition of a determinant, that 

det(mP,Q)=mdet(P,Q)=det(P,mQ) (I.) 

Suppose now that P and Q after an interval of time At become respectively (P)-f-(AP), 
(Q)+(AQ), the sign + here denoting the complete sum. Then the complete increment 
of det (P, Q) is 

det (P+ AP, Q+ AQ)-det (P, Q)= l 

det (P+ AP, Q+ AQ)-det (P+ AP, Q)+det (P+ AP, Q)-det (P, Q).J ( L ' 

But it follows from formula (II.) of section 26, that 

det(P+AP, Q+AQ)-det(P+AP,Q)=det(P+AP, AQ); 
and similarly, 

det (P+ AP, Q)-det (P, Q)=det (AP, Q). 
Substitute, then, these values in the above equation (II.), and divide both sides by A^by 
means of the above formula (I.), and finally let At diminish without limit. We thus 
obtain for the complete differential coefficient of det (P, Q) 

det (P, T> t (Q))+det (D,(P), Q). 
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We have therefore the equation 

D,{det (P, Q)} =det {P, T> t (Q)} +det {D,(P), Q[ . 

The same equation may also be proved by considering the algebraical determinants 
which represent the projections of det (P, Q) ; and it may in fact be easily deduced from 
the following identical equation, 

dtW'Pi 9*P*)—\dt P* dt P') + \S* dt Sy dt) 

29. It may be here observed that the formula? (I.) and (II.) in section 26, and the 
formula? (I.) and (II.) in section 28, show that there exists an intimate symbolical connexion 
between det (P, Q) and the product P, Q. In fact the only difference between their 
symbolical properties consists in P and Q not being commutative in the expression 
det (P, Q), and being so in the expression for the product. 

30. There is one more proposition which is often very useful in analytical dynamics. 
Let it be required to find det (R, Q'), where Q' itself equals det (P, Q). Let the required 

line det (R, Q') be denoted by U. Then, by the definition of a determinant, U is perpen- 
dicular on R and on Q', which last line is itself perpendicular on the plane containing 
P and Q. Hence it follows that U is perpendicular on R and in the plane containing 
P and Q. 

We have still to find the magnitude of U. For this purpose let the angle which R 
makes with the plane containing P and Q be \|/, so that \f/ is the complement of the 
angle between R and Q'. 

Moreover, let 6 be the angle between P and Q, and let the magnitudes of P, Q, Q', R be 
denoted by^?, q,q', r respectively. Then, since U is det (R, Q'), it follows from the defini- 
tion of a determinant, that the length of U equals rq'sia (o - ^) or r S' cos ty Simi- 
larly, g l =pq sin 6. Hence the length of U equals pqr sin 6 cos \J/. 

There are two cases especially which frequently occur in dynamics, first, when R is 
identical with Q, and secondly, when R is perpendicular on Q. 

Let us first take the case of R being identical with Q ; then \J/=0 and r=q. There- 
fore the required determinant is a line in the plane containing P and Q, and perpendi- 
cular on R or Q, and its length equals pg 2 sin 6. 

If, moreover, Q is perpendicular on P, then the required determinant is in the 

direction of P, and its length equals pq 2 , since d=z~. So that, if P is perpendicular on 

Q, we see that det {Q, det (P, Q)} is a line^<f in the direction of P ; and therefore evi- 
dently det {Q, det (Q, P)} is a linej^ 2 opposite to P. 

Let us now take the case of R being perpendicular on Q. Then it might be easily 
proved by spherical trigonometry, that sin 6 cos \|/ equals the cosine of the angle between 
R and P. But we will prove this by analysis, because in doing so we shall meet with 
formulse which will be of use in the sequel. 



482 ME. A. COHEN ON THE DIEEEEENTIAL COEFFICIENTS 

Let, then, the components of P, Q, Q', E parallel to any three axes of coordinates 
be denoted by p x ,p y ,p z , q x , Sec, q x , Sec, r„ Sec Then, if we denote the components of 
U=det(E, Q') by u x , u y , u x , we have, by section 24, 

w x =q z r y —q y r z ; • (I.) 

and since Q'=det (P, Q), we have 

q y =q x p-qg x , 

Hence, substituting the values of q z and q' y in (I.), we get 

^■=^{.q y r y +q^ z )—q x {p y r y +p z r z ). 

Now by hypothesis E is perpendicular on Q ; hence 

q,r x +q y r y +q z r z =0 ; 
therefore 

XT V "T~ jt z z^~ """ q%' %' 

Therefore 

v*= — qx(p*r x +p y r y +p z r z ). 

But as q and r denote the magnitudes of P and E, it is evident that 

Pfm +Py^y +iV.=3» , COS <p, 

where <p denotes the angle between P and E. Therefore 

u x =— prq x cos <p. 

Similarly 

u y =— prq y cos <p, 
u z = —prq z cos <p. 

Therefore the line U of which u x , u y , u z are the components is a line in direction oppo- 
site to Q, and whose length equals prq cos <p, q being the length of Q. Hence if E be 
perpendicular on Q, then det (E, det (P, Q)) equals — pqr cos <p in the direction of Q. 

It follows from the above proposition, that, if Q' or det (P, Q) represent a force or 
acceleration which acts at the extremity of the radius vector E, and if Q be perpendicular 
on E, then the moment-axis of that force or acceleration about the origin is —pqr cos p 
in the direction of Q, and the moments of such force or acceleration about the coordi- 
nate axes are respectively — pr cos <p q x , —jpr cos <p q y , —pr cos <p q z ; and as <p is the angle 
between P and the radius vector E, pr cos <p=xp x -\-yp y -\-z]) z , if x, y, z be the coordinates 
of the extremity of the radius vector. 

Chapter III. 

31. We are now in a condition to treat fully of the motion of a particle in space of 
three dimensions ; and it will be found that the propositions which have just been 
proved concerning the determinants of lines, will enable us to show how all the results 
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arrived at as to a particle's motion in one plane may be extended to motion in space 
generally. 

32. Suppose Q to represent a line drawn from the origin, varying both in direction 
and magnitude in any manner whatsoever, and let it be required to investigate D^Q) the 
complete differential coefficient of Q. 

Let the length of Q be q at time t, and let the direction of Q be revolving at time t 
about a line whose direction is that of the line represented by O, and let the length of 
Q, be the angular velocity nr, with which Q's direction is revolving at time t. 

It has been already shown in section 12 that D^Q) is in all cases the complete sum 
of the two partial differential coefficients which are obtained by varying separately the 
length and direction of Q. Now the former partial differential coefficient is evidently 

Iji II to Q], and the other partial differential coefficient is, by section 25, equal to 

det (Q, Q). Hence we have the following fundamental equation, 

D,(Q)=(§||toQ)+det(Q,Q) (I.) 

33. It is not difficult to deduce from the last equation the expression for D|(Q), the 
second complete differential coefficient of Q. In order to find that expression we must 
take the complete differential coefficient of each of the expressions of which the right- 
hand member of equation (I.) is composed. For this purpose represent for a moment 

the line (■£ [| to Q) by Qj. Then it follows from the fundamental formula of the pre- 
ceding section, that 

Df(Q)=D,(Q 1 )+D, (det (Q, Q)). 

Now the formula (I.) of the last section gives evidently 



D ( (a)=(^ II to Q,) +det (Q, Q,), 



or 



(JlJtoQ^+det^QJ. 

Moreover we have, according to section 28 of the preceding Chapter, 

D,{det (Q, Q)} =det {D,(Q), Q} +det {&, T> t (Q)\ . 
But since 

D,(Q)=(Q 1 )+det(Q,Q), 

it follows from section 26 of the preceding Chapter, that 

det (Q, D,(Q))=det (Q, Q,)+det {Q, det (Q, Q)} . 

Therefore, collecting the above results, we obtain 

D?(Q)=D i (Q 1 )+D < {det (Q, Q)} = 

( % || to q) +2 det (Q, QJ+det (D,(Q), Q)+det^Q, det (Q, Q)[ . 
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The two last terms of this expression are evidently what would be D*(Q) if 
^ were zero, that is to say, if Q did not vary in magnitude; and ( -p || to Q J is evi- 
dently what D|(Q) would be if Q did not vary in direction; so that we have the 
following proposition : — 

" D|(Q) is the complete sum of the two partial second differential coefficients obtained 
by varying separately the length and the direction of Q, together with 2 det (Q, Q x ), 

where Q, is the line (^ \\ to <A" 

34. Suppose Q to be E the radius vector of a moving particle, the length of which 

radius vector is r, then D?(Q) is the particle's acceleration; Q : is l-n || toE), and 

is therefore the velocity along the radius vector. If, then, we denote this by E 15 the 
equation arrived at in the last section shows that the acceleration is compounded of 

^||toE)+2det(Q,E 1 ), 



and of what would be the particle's acceleration if E did not vary in magnitude, that is 
to say, if the particle simply revolved about the origin. And this latter acceleration is 
again compounded of det (D ( (Q), E) and det {Q, det (Q, E) j- . The last line is, by section 
30, a line drawn from the extremity of E, or from the particle, perpendicular to and 
towards Q, and whose magnitude is ra- 2 ^, p being the length of that perpendicular. 

35. The above result is, however, but a particular instance of the theory of the motion 
of a particle relatively to axes which revolve about the origin, a subject which we are now 
in a condition to treat of very simply in its utmost generality. That theory will be found 
to depend upon the solution of the following problem : — 

" Supposing the axes of coordinates Ox, Oy, Oz to revolve round the origin O about 
an axis Q at time t with angular velocity is (which is the length of Q), it is required to 
find the complete differential coefficient of a line Q, the components of Q along the 
coordinate axes being given." 

Let Q have for its components Q x , Q y , Q s , and let the respective lengths of these be 
q x , q s , q z . Then evidently 

Q=(Q.)+(Q r )+(Q»). 
Therefore 

But, by the fundamental formula of section 32, we have 



D*(Q*)= (% II to 0*) +det (Q, QJ, 

D ^)=(^il to0 ^)+ det ( Q '^)' 
D,(Q,)= (§? || to Oz) +det (O, a>' 
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Now we know, from section 26 of Chapter II. , that 

det (Q, QJ+det (Q, Q y )+det (12, Q z ) 

=det{Q,(Q,)+(Q,)+(Q,)* 
=det (Q, Q). 
Therefore 

D,(Q)= (§= i! to Oar) + (I 1 II to Oy) + (§ || to 0«) +det (Q, Q). 

This formula is true whether the axes be rectangular or oblique, and may be made 
the basis of all the formulae of relative motion. 

It may be observed that, if the coordinate axes did not move, D,(Q) would be equiva- 
lent to 

(^ || toO.) + (f |, to Oy) + (# || toO.) . 

So that the line represented by the last expression may be called the differential coeffi- 
cient of Q relatively to the moving axis, or, more briefly, the relative differential coefficient 
of Q. The above formula of the last section therefore shows that the complete differential 
coefficient of Q is the relative differential coefficient of Q together with det (Q, Q). 
This proposition exactly corresponds with the proposition in section 17 of Chapter I. 

36. Assuming now the axes of coordinates to be rectangular, we know, from formulae 
(II.) in section 24 of Chapter II., that det(Q, Q) has for its components 

g^y—ciy^z parallel to Ox, 
q x vT z —q z vT x parallel to Oy, 
q 1/ vT x —q x vT !/ parallel to Oz. 

Therefore it follows from the preceding section, that the components of T> t (Q) are 

<% . 
dq* 

These formulae are in fact simply the analytical expression of the fundamental propo- 
sition in the preceding section, and correspond exactly to the formulae in section 18 of 
Chapter I. 

37. Let us now apply the above formula? to dynamics, and first to the velocity of a 
particle. 

Suppose, then, a particle to move in space in any manner whatsoever, and suppose that 
the rectangular axes of coordinates revolve about a line Q at time t with angular velocity 
nr, zs being the length of Q. Let v x , v y , v z be the components of the particle's velocity, 
and zs x , ©-J,, vr z the components of Q. Then, since the velocity is the complete differential 
coefficient of the radius vector K, of the particle, and since x, y, z are the components of 

mdccclxii. 3 u 
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E, it follows at once from the formulae of the preceding section, that 

dos 

V v=di+ X,BT *— Zv >z» 
dz , 

These formulae simply express the fact, that the absolute velocity is equivalent to the 
relative velocity together with det (Q, K). 

38. Let us next apply the formulae to the acceleration of a particle. Let, as before, 
v x , v r v z be the components of the particle's velocity V, and let f^f^f,. be the components 
of the particle's acceleration. Then, since the acceleration is the complete differential 
coefficient of the velocity, of which v x , v y , v z are the components, it follows at once from 
section 36, that 

n dv x 

Jx—-dt+ V ^y — V y ™z, 

j, dv z 

If we substitute in the last equations the values obtained for v x , v p v z in the preceding 
section, we obtain 

j, (Px dar„ «fer r , „ /dz dy \ , , . , , 

/«=#« +z ~~dr~y~dT+ 2 '\Tt r *y~Tt' l!r *) +{yw«— ««■»)'»,— (Xm M —Z*r.) m n 

and similar formulae for^ and/^. 

These are the ordinary formulae. It would not be difficult to deduce their real 
meaning from their analytical form ; but it will be better first to prove the result of 
such interpretation in a different and more direct manner. 

39. We have already seen that, if V denote the particle's absolute velocity, and R the 

dy (Ln dz 

radius vector, and if Y 1 denote the relative velocity which has -77, ■£, ~r, for its compo- 
nents, then 

V=V 1 +det(Q,E). 

Let then F denote the particle's absolute acceleration, and let Y^ denote the particle's 

vL of* d 11 (l z 

relative acceleration which has -p-, -p-, -p- for its components ; then 

F=D < (V)=D < (V 1 )+D,det(Q,R) (I.) 

Now it follows from the fundamental proposition in section 35, that 

D,(V I )=F 1 +det(Q,V 1 ). 
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Moreover we know from section 28, that 

J) t det (Q, R)=det (D,(G), R)+det (Q, D,(R)> 

But since D,(R)=V = V 1 +det (O, R), therefore 

det (Q, D,(R))=det (Q, VJ+det {Q, det (Q, B)}. 

Hence, collecting the above results, and substituting them in the equation (I.), we 
find 

F=F 1 +2 det (Q, V,)+det (D,(Q), R)+det{Q, det(Q, R)f . 

It may be observed as to this formula, that if V, = 0, that is to say, if the particle had 
no relative motion, and moved as if rigidly connected with the axes of coordinates, then 
the two first terms of the last equation would vanish ; and therefore its other two terms 
are what the acceleration would be if the particle had no relative motion, and they 
represent what may therefore be conveniently termed the particle's system-acceleration. 
French writers have given to this acceleration the name of "acceleration d'entrainement; " 
it is the acceleration of a point which is in the position of the moving particle, and 
which is supposed to be rigidly connected with the system of moving axes, and I there- 
fore propose to call it " system-acceleration." Using then this expression, we have the 
following proposition : — " The acceleration of a particle is equivalent to its acceleration 
relatively to a system of axes revolving about a fixed point, together with the system- 
acceleration corresponding to the particle and together with an acceleration equal to 
2 det (Q, V,), V! being the particle's relative velocity, and Q the axis about which the 
system is revolving at time t." Or, more briefly, a particle's absolute acceleration equals 
the complete sum of its relative acceleration and of its system-acceleration together with 
2 det (Q, V,). Such is the brief expression of Coeioli's beautiful and very useful pro- 
position concerning relative motion. 

40. We have just seen that the particle's system-acceleration is compounded of 
det (D t (Q,), R) and det \Q, det (Q, R)}. As regards the latter line, it is clear from 
section 30 that it is in the direction of the line drawn from the particle perpendicular 
to and towards the axis Q, and that its magnitude is v^p, p being the length of that per- 
pendicular. It is therefore equal and opposite to what is usually called " the centrifugal 
force." 

As regards the other line det (D t (Q), R), it is of course at once determined as soon as 
D 4 (Q) is known. Now if nr x , ar y , m s be the components of Q, it follows from the funda- 
mental proposition in section 35 that D f (Q) is equivalent to 

(** 1| to 0*) + (^- 1| to Oy) + (*£ || to O*) +det (Q, Q). 
But it is clear, from the very definition of a determinant, that det (Q, Q) is zero.: 
Hence we see that the components of D f (Q) are -~, -^, -~. This is an important pro- 
position often used in the dynamics of a rigid body, and generally proved by means of 
a good deal of analytical work. It is usually expressed in the following manner. If 

3u2 
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w„ rar 2 , vr 3 be the components of w along fixed axes, and ts m vr y , gt, be its components along 

moving axes which coincide with the former at time t, then J^l = ^, -^=J^ Jp = _?5. 
& ' dt. dt dt dt' dt dt 

It is evident that this amounts to saying that D,(Q) has ~, ~ y , -^ for its components ; 

Ctl G/l CtZ 

and we have just seen how that proposition follows at once from the fundamental 
theorem in section 35, and from the self-evident fact that D,(Q, Q)=0. 

41. Recapitulating then the results of the two last sections, we see that a particle's 
system-acceleration is equivalent to det (D^(Q), R) minus the centrifugal force, and that 
the absolute acceleration of the particle is compounded of the relative acceleration of 
the particle, its system-acceleration, and 2 det (0, VJ, Y Y being the particle's relative 
velocity. 

If we now look back on the analytical expressions obtained in section 38 for the com- 
ponents of the absolute acceleration, it will be easy to see their full meaning. The 

expression -p is the component of the relative acceleration. The expression z-^ — y -p 
is the component of det (D,(Q), R), since, as we have seen, 10,(12) has for its components 
-^, -p!, -^. The expression 2 I Tt' m y~dt' me ) * s ^ ne com P onerL t of 2 det (Q, VJ, since V! 
has for its components -^-*, -?pf, -p*. Finally, it may be easily shown by analytical geo- 

OjZ tit ttZ 

metry, that the expression 

(&*** — ^yl^y ~ ( XvT * — z ™x)™z 

is the component of the line is*p drawn from the point (x, y, z) on the line whose direc- 
tion-cosines are proportional to m x , w y , -m„ p being the length of that perpendicular. 
Hence it is manifest that the analytical formulae in question merely express the proposi- 
tion enunciated at the commencement of this section. 

It may, finally, be observed that the above results might also have been easily deduced 
from the formula in section 34 for the acceleration along the radius vector in exactly the 
same manner as the corresponding analytical formulae for the relative motion of a par- 
ticle in one plane were deduced in section 22 from the formula for the acceleration of 
the particle along the radius vector. 

42. If the origin of coordinates also moves, it is evident that the particle's actual 
acceleration is the resultant of the acceleration of the origin and the acceleration rela- 
tively to the origin. Hence substituting for the latter acceleration the expression 
already found for it, it is easy to see that the particle's actual acceleration is, as before, 
the resultant of the relative acceleration, an acceleration represented by 2 det (O, VJ, 
and the particle's system-acceleration, but that the system-acceleration is now the 
resultant of the acceleration of the origin, and of the system-acceleration relatively to 
the origin, for which latter system-acceleration we have already obtained the expression. 
Now in whatever way a system moves, the motion may be decomposed into a motion of 
translation and a motion of rotation. Hence we see that a particle's absolute accelera- 
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tion is in all cases the resultant of the relative acceleration, the system-acceleration, and 
an acceleration equal to 2 det (Q, V), where Q is the axis about which the system is 
turning at the time, and V is the relative velocity of the particle. 
This is the most general form of Corioli's theorem. 

43. One of the most important illustrations of the theory of relative motion is the 
motion of a heavy particle relatively to a system which revolves uniformly about a fixed 
axis ; for this includes the case of a falling body and the pendulum, where the earth's 
motion is taken into account. 

Suppose then a particle of mass m to have for its actual weight W, and for its appa- 
rent weight W, so that a force — W would keep the particle in relative equilibrium or 

/W'\ /w\ : , 

apparently at rest. Then evidently ( — J — - 1 — 1 is equivalent to the particle's system- 
acceleration. 

Let then the particle be acted on by a force P over and above the weight W, and let 
the particle's actual acceleration beF, its relative acceleration F 15 its system-acceleration 
F 2 . Then clearly 

*=(?) + (5\ 

\mj ■ \rn J 

But by Corioli's theorem 

F=(F 1 )+(F 2 )+2det(Q,V); 

/W'\ /W\ 
and we have just seen that the system-acceleration F 2 = ( — j — I — ] . Hence it follows 

that 

(F,)+2det(Q,V)=@ + (^): (I.) 

W . 

— is the apparent acceleration of gravity, and is generally denoted by g. 

44. The above formula is quite general ; but in most cases g may be considered as 
constant both in magnitude and direction, its direction being the vertical direction at the 
point of reference or origin. 

We have then the formula 

F i=®+(#)- 2det ( Q > V ') W 

This simple formula enables us to solve easily all problems concerning the motion of 
a heavy particle relatively to a spectator on the earth. The formula shows that the 
relative acceleration is found, just as if the earth did not move, by substituting the 
apparent for the actual force of gravity, and by adding on a force —2m det (Q, V,), 
where Y x is the particle's apparent velocity. 

45. Let us take the vertical downwards as axis of z; let the axis of xhe the horizontal 
line drawn from north to south, and let the axis of y be the horizontal line drawn from 

west to east. Then the equation to the earth's axis is evidently 7^r^ = ^"^5 ^ ^ denote 
the latitude of the spectator's position. 

Therefore «r,=w cos X, w„=0, ^rrro-sinX. 
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Moreover we know that the components of det (Q, V,) are 



dz dy . dy 

doc dz (dx . dz 

di™* 



dz (dx . dz \ 

-Jt^=^{-dt smX -Jt COsK )- 



dy dx dy 

Let then the force P, which, besides gravity, acts on the particle, have for its 
components X, Y, Z, then, as F„ the relative acceleration, has for its components 

<jp x ,jpy $l z ; 

^2, ^p, -^2, it evidently follows from equation (II.) of the preceding section, by revolving 

along the axes, that 

d z x X dy 

dT* = m + 2 ™ SlnK W 

d*y Y , „ (dz dx . \ /TTT . , 

A» = ™ + 2w (*. C °8X-^SUlXJ, V ....... (HI.) 

^ Z , _ dy 

ffl* fllJ fl'Z 

On multiplying these equations by ~rr, -±, t- respectively, and atlding, we find 

dx d?x dy d^y dz <Pz X dx Ydy Zdz dz 

di ~dP>di df~>~Tt dF = m di^"m Tt'm It'^ ~dt 

This equation may also at once be deduced from the formula (II.) if we resolve 
along the direction of the particle's relative motion, and observe that det (Q, V,) is per- 
pendicular on that direction which coincides with the direction of V,. 

By integrating the last equation, we see that the equation of vis viva applies to the 
particle's relative motion just as if the particle's relative motion were its actual motion, 
with this difference only, that for the actual force of gravity the apparent force of gravity 
must be substituted. 

If the particle be a free particle acted on by no forces but gravity, then X=0, Y=0, 
Z=0, and the equations (III.) are linear, and are therefore easily integrated. 

Moreover if v, be the relative velocity, the equation of vis viva gives 

v 2 1 =2g(2— h), h being a constant. 
46. If the particle be suspended by a string from a point fixed to the earth, then if 
that point be taken as the origin, and X, Y, Z be the components of the string's tension, 
we evidently have Zy— Yz=0, X.z—Zx=0, Yx— X^=0; and substituting those values in 
the preceding equations (III.), we shall obtain two independent equations, which, together 
with the equation x 2 -\-y*-\-z 2 =& constant, will determine the particle's relative motion. 
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But those resulting equations can be found far more simply and directly in the following 
manner. 

For this purpose let us revert to the fundamental formula 



F 



'=@+(£)- 2det (0>V) (I.) 

Now section 30 of Chapter II. shows us how to find the moment-axis with respect to 
the origin of det (12, V). If, namely, ra- and v be the magnitudes of 12 and V, and p be 
the angle between 12 and the radius vector, then the moment-axis of det (12, V) is 

—mvr cos p. Therefore the moment about the axis of x of — 2 det (12, V) is 2srr cos p -^; 

dij 
and similarly, its moments about the axis of y and z are respectively 2wr cos p -j t and 

Hz 
1-mr cos p j.. Moreover, since p is the angle between the radius vector and 12, it is evident 

that 

vxr cos p=vrjc-\-v? !/ y-\-vr !: z=:x>s(x cos A+iZ sin X), 

X being the latitude of the origin. 

Let us now take the moments of Fj about the axes of coordinates, and equate them to 
the moments of those components of F x which are given in the formula (I.). 

P 

The moments of - about the axes are zero in this case of the pendulum. The 

til 

moments of g about the axes of x, y, and z are respectively gy, —gx, and ; and the 
moments of —2 det (12, V) we have just found. Hence we obtain at once the following 
three equations : — 



d z z d*y . „ dx . . . 

y-d?- z ~di= ^+2*r^(#eosX+zsmA), 

z ■^—x-^=—gx+2nr J- (X cos X-f 2 sin X), 
d?y cPw n dz , . . , 

cc 'df-yi? = 2w Tt ( x cos x+ z sm x )- 



(II.) 



These are the three equations given in Hansen's elaborate ' Theorie der Pendel- 
Bewegung,' and which are generally obtained by means of very complicated analysis. 

One simple equation can be deduced by means of the proposition contained in section 
45 ; for the principle of vis viva gives 



;sy+ ($y+(S) ■=*.)+.«—«. 



Moreover af+y^+z* is a constant, and these two equations, combined with any one of 
the equations (II.), determine the particle's motion. 
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Chapter IV. 

47. As soon as we pass from the statics or dynamics of a particle to the statics or 
dynamics of a system of particles or of a rigid body, we find that two forces which are 
equal and parallel to one another are not equivalent to one another, and that we have 
to take into account the position as well as the magnitude and direction of a force. 
Notwithstanding this, we are enabled by means of an elementary principle of statics to 
confine our operation and our notation to lines passing through one and the same point. 
For suppose a force P to act at a point to of a rigid body, and apply at the origin O two 
equal and opposite forces P and — P, then P at m is equivalent to P at O and the couple 
whose forces are P at to and — P at O. Let the axis of this couple be denoted by G; 
the couple, being completely determined by G, may be called the couple G. It is 
extremely convenient to have a name for the line G, indicating briefly its connexion with 
the force P at to, and I shall adopt that given to it by French writers*, and shall call G 
the moment-axis about O of the force P at m. 

It has been proved in section 25 of Chapter II., that the line G, being the axis of the 
couple whose forces are — P at O and P at to, is equal to det(— P, R), where E is the 
radius vector of the particle. Hence we have 

G=det (-P, R)=det (R, P). 
We thus see that the force P at to is completely represented and determined by the 
two lines P and G drawn from the origin, G being the moment-axis with respect to the 
origin of P at to, and being equal to det (R, P). 

48. Suppose now that we have a system of forces P„ P 2 , &c. acting respectively at 
points TOj, to 2 , &c. of a rigid body. Then it is clear from statics that the given system 
of forces is equivalent to a force P at the origin O and a couple whose axis is G, where 
P is the complete sum of the forces P„ P 2 , &c. supposed to be collected at the origin, 
and G is the complete sum of the moment-awes (G„ G 2 , &c.) (about the origin) of the 
forces Pi at to,, P 2 at to 2 , &c. 

49. We will now apply the above considerations to dynamics. Since the acceleration 
is the complete differential coefficient of the velocity, it is evident that the line which 
represents the moving force of the particle is the complete differential coefficient of the 
line which represents the particle's momentum; or, more briefly, the moving force is the 
complete differential coefficient of the momentum. 

Let now P represent the moving force, and U the momentum of a particle to, P and 
U denoting straight lines ; then, if we treat the moving force and momentum as if they 
were statical forces, it is clear that P at to is equivalent to P at the origin O and a 
couple G, where G is the moment-axis about O of P at to ; and similarly, the momentum 
U at to is equivalent to U at O and a couple of momenta whose axis is H, where H is 
the moment-axis about O of U at to. We have just seen that P is the complete differ- 
ential coefficient of U, and we will now prove that in like manner G is the complete 

* See Delattnay's ' Mechanics,' page 254. 
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differential coefficient of H. If, namely, R denotes the radius vector of the particle, 
G=det (E, P), and similarly H=det (R, U). Now, if we differentiate the last equation, 
H=det (R, U), we obtain, according to section 28, 

D,(H)=det(K,D,(U))+det(A(R),U) (I.) 

But D f (R) is identical with the particle's velocity, and is therefore in the direction of 
the momentum U. Whence it follows, from the very definition of a determinant, that 

det (D t (R), U)=0. 

Therefore the above equation (I.) becomes, since D,(U)=P, 

D,(H)=det (R, D,(U))=det (R, P)=G. 

This is an important result. It shows that the moment-axis about any point of the 
moving force of a particle is the complete differential coefficient of the momentum, and 
that therefore the moment of the moving force about any line is the differential coefficient 
of the moment of the momentum. 

The above result may also be easily deduced from the identical equation 



d z y d*x df dy dx\ 



50. The proposition which we have just proved may be easily extended to a system 
of moving forces and momenta of the particles of a rigid body. For, according to section 
43, the system of moving forces is reducible to a moving force at the origin O, and a 
couple G. And the system of momenta may be similarly reduced to a momentum U, 
and a momentum-couple whose axis is H. Now we have seen that P is the complete 
sum of the moving forces, and that each moving force is the complete differential coeffi- 
cient of the corresponding momentum. It therefore evidently follows that P is the 
complete differential coefficient of the complete sum of the momenta, or of U. Hence 
P=D,(U). Moreover we have seen that G is the complete sum of the moment- 
axes about O of the moving forces, and that each of these moment-axes is the complete 
differential coefficient of the moment-axis of the corresponding momentum. Hence it 
follows that G is the complete differential coefficient of the complete sum of the moment- 
axes of the momenta. Hence G=H.(H). 

This result may be also easily proved by means of the identical equations 



_, / cl?x\ d _ / dx\ 

_, / d?y d?x\ d _, / dy dw\ 



51. The science of the dynamics of a rigid body is founded upon D'Alembert's prin- 
ciple, which asserts that the moving forces of a body's particles are together statically 
equivalent to the impressed forces acting on the body. If therefore these external forces 
be reduced to a force P at the origin O and a couple G, then P and G are equal 

mdccclxii. 3 x 
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respectively to what was denoted in the preceding section by P and G ; and we have 

therefore 

P=D ( (U), G=D,(H). 

In other words, if we treat the momenta as statical forces, and reduce the system of 

momenta of a body's particles to a momentum U at O and a momentum-couple whose 

axis is G, then the external forces acting on the body are equivalent to the force D^U) 

at O and the couple of forces whose axis is D,(H). 

Since G=D/H), the resolved part of G along any fixed line will be the differential 
coefficient of the resolved part of H along that line ; or, in other words, the sum of the 
moments of the external forces about any line equals the differential coefficient of the 
sum of the moments of the momenta of the body's particles. 

The above results may be easily deduced from the ordinary equations 

S(X)=s(mg)=^(m|),&c, 

But it will be generally found far better not to use those six equations at all, and 
simply to bear in mind the fact which they alone express, namely, that P =D t (U), 
G=D,(H). 

52. It will be convenient to recapitulate once for all the notation and phraseology I 
shall constantly use in the sequel. The system of momenta of a body's particles, or 
what may be called the body's momenta-system, is reducible, if we treat the momenta as 
forces, to a momentum at a point O, and a couple of momenta. The former I call the 
body's momentum, and denote it by U ; the latter I call the body's momentum-couple about 
O, and denote its axis by H. U and H may both be represented by straight lines 
through the origin O. It is to be observed that U remains the same wherever O be 
taken, but that H changes with the position of O. 

The components of U and H along the axes of coordinates will be denoted by 
Ua,, JJ r U,,, H a , ~H. y , Hs respectively, the magnitudes of all these quantities being repre- 
sented by the corresponding small letters. Thus h x will be equal to the resolved part of 
H along the axis Ox, and will therefore equal the sum of the moments of the momenta 
about Ox. 

The system of moving forces is reducible to a moving force D,(U) at O, and a couple 
whose axis is D^H) ; and it follows from D'Alembert's principle that, if the external 
forces be reduced to a force P at O and a couple G, then 

P=D ( (TJ), G=D,(H). 

53. The next step will be to investigate the expressions for U and H. 

In the first place, U can be easily found. For, let R denote the radius vector of a par- 
ticle of mass m; then, if 2 denote the operation of taking the complete sum of lines, we have 

U=2mD,(E) 
=D£(mR). 
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Now, by a well-known proposition, it is evident that 2j(mB,)=MR, where M is the mass 
of the body and R is the radius of the body's centre of gravity. Hence 

U=ML\(R); 

and therefore, if V denote the velocity of the centre of gravity in magnitude and direction, 
U=MV, or, in other words, U, the body's momentum, is the momentum of the bodys mass 
collected at the bodys centre of gravity. 

54. In the next place we have to find H. The investigation will be very much 
facilitated by the following consideration. A body's motion is said to be compounded 
of motions «, j3, y, if the velocity of each of the body's particles may be considered as 
the resultant of the respective velocities due to the motions of a, f3, y separately. In 
such case the momentum of a particle will evidently be the resultant of the momenta 
due to each of the motions «, (3, y separately ; and since the resultant of the momenta 
of all the particles will be the same in whatever way we group them together, it is 
evident that we have the following proposition : — 

" The resultant of the momenta of a body's particles, or the body's momenta-system, is 
the resultant of the momenta-systems due to each of the motions «, /3, y." 

Thus a body's motion may be decomposed into a motion of rotation and translation. 
Hence the body's momenta-system may be found by compounding the momenta-system 
due to the motion of rotation with that due to the motion of translation. 

Again, a motion of rotation may be decomposed into rotations about three axes. 
Hence the momenta-system of a body which rotates about a fixed point is the resultant 
of the momenta-systems respectively due to the separate motions of rotation about the 
three axes. 

55. Let us then first investigate H for a body having simply a motion of translation. 
Let v be the velocity of translation in the direction of a line AB at time t, then the 

momentum of a particle of mass m is mv in the direction of AB. Hence the momenta- 
system consists of a number of momenta parallel to one another, and proportional to the 
masses of the respective particles. Their resultant is therefore Mv at the centre of 
gravity, M being the body's mass. Hence the momenta-system is reducible to Mv at 
the centre of gravity. Therefore H, the axis of the body's momentum-couple about O, 
is the moment-axis about O of Mv at the centre of gravity, and is zero, if the point O 
coincides with the centre of gravity. In the latter case, since H=0, therefore D<(H) = 0, 
therefore G=0, or the moment of the external forces about any line through the centre 
of gravity is zero for a body which has simply a motion of translation. 

56. The next simplest case is that of a body rotating about a line. Take that line as 
axis of z, and suppose the body of mass M to be revolving at time t about that line with 
angular velocity «-.. It may be easily shown in the ordinary way, that the sum of the 
moments of the momenta about the axes of x, y, z are respectively 

— srSfwm), —Gf%(myz), cr2(mr 2 ). 
3x2 
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So that, using the notation of section 52, we have 

h x =— 7s%(mxz), \—— vr%(myz), h g = ^(mr*). 

If the axis of z be a principal axis, we have 

%{mxz) = 0, t(myz) = 0, 
therefore 

4=0, \=0, 4=w2(mr 2 ); 

and consequently H is a line in the direction of Oz, and equal to the product of sr and 
the moment of inertia about Oz. 

57. If the axis about which the body rotates is neither a fixed line nor a principal 
axis, it is more convenient to express II, the body's momentum-couple, in the following 
manner. 

The body's rotation about the point O may be considered as compounded of motions 
of rotation about the principal axes at O. Let «r„ w 2 , m 3 be the angular velocities of those 
component rotations, and let A, B, C be the respective moments of inertia about the 
principal axes. We have shown in the preceding section that the body's momentum- 
couples due to the three separate rotations about the principal axes would have for their 
respective axes lines along the principal axes and equal to Ato-j, Bw 2 , CV,. It follows, 
therefore, from section 54, that H, the body's momentum-couple, is the resultant of the 
three couples, whose axes are repectively Ato-,, IW 2 , CW 3 . In other words, H, the axis of 
the body's momentum-couple, has Asr,, IW 2 , (V, for its components along the principal 
axes. 

58. The results of the last two sections may be also proved in the following more 
direct manner. 

Take any rectangular axes as axes of coordinates. Let ■&„,, ■m p vr z be the components, 
along those axes, of the body's angular velocity of rotation. Let v x , v s , v z be the compo- 
nents of the velocity of a particle of mass m, whose coordinates are x, y, z. 

If then h m h r h z be the components of H, we have 

h x =tm(yv z —zv y ). 

But V z =yi<r x — XiXy, « y =^w s — Zta x . 

Therefore h x = ■as£m{y l + z 2 ) — ■n y %{myx) — rs£(mxz) . 

If, then, the moments of inertia about the axes of x, y, z be denoted by A, B, C respect- 
ively, and if we denote ~2(myz) by A', 1{mxz) by B', %(myx) by C, the last equation 
becomes 

K— Aw,— CW,— BW, ; 

and similarly, ^B^-AWj-CW,, 

\ = CVj — BV, — AWy 

These results are true, whatever rectangular axes of coordinates be taken ; but if they 
be principal axes, then A'=0, B'=0, C'=0, and therefore we have, as before, 

h x =Anr x , A,=Bsr y , h z =Cnr z . 
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59. We are now in a condition to solve easily the problem of the motion of a body 
rotating about a fixed line, or about a fixed point under the action of any forces. 

First, let us take the case of a body of mass M revolving about a fixed line. Take 
that line as axis of z, and choose for the axes of x and y any lines fixed in the body 
which are perpendicular to one another and to Oz. 

It was proved in section 52 that H, the axis of the body's momentum-couple, has for 

its components 

h x = — GT%(mxz), h y =—m%(myz), h z =M.lc i tjj, (I.) 

M# 2 standing for the moment of inertia about Oz. 

Moreover we know from section 48 that U, the body's momentum, is MV, where V 
is the velocity of the centre of gravity. Now if x, y, z be the coordinates of the centre 
of gravity, V has evidently for its projections on Ox, Oy, Oz, — my, izx, respectively. 
Hence the components of U are equal to 

u x =— M.vjy, Uy=M.mx, u e =Q ; (II.) 

knowing, then, the components of U and H, we can easily find the components of D^U) 
and L\(H). Using the notation of section 52, P=D < (U) and G=D,(H), and the com- 
ponents of P and G may be denoted by ~P X , P r P 2 , and G„ G y , G z respectively. 
In the problem now before us, the axis of z does not move ; hence evidently 

P.=|(«.) and G~(h.) (III.) 

But as to the axes of x and y, they revolve about the axis of z with an angular velocity 
■m at time t. Hence by the elementary formulae of section 18 in Chapter I., we have 

Similarly, i (IV.) 

Ga,=J t (h)— hj/nr, G s =j t (\)-\-h x iB. 

Let, then, external forces acting on the body be reduced to a force at O whose compo- 
nents are X, Y, Z, and to a couple the components of whose axis are L, M, N. Let 
the reactions of the fixed axis be similarly reduced to a force whose components are 
X', Y', Z', and to a couple the components of whose axis are L', M'. Then, by D'Alem- 

beet's principle, 

X+X'=P„ &c., 

L + I/=G„&& 

Therefore, substituting the values (I.) and (II.) in equations (III.) and (IV.), we obtain 
the following six equations : — 

X+X'= -My ~-Mv% 

Y+Y'= MS^-Mw^, 
Z+Z'= 0, 
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L+L' = _ t(mxz) -j t J r zs i t(myz), 
M+M'= — 2(myz) -^ —vr 2 2(mxz), 

These six equations are those ordinarily given in text-books, and their full import and 
meaning is now apparent. The first three have for their right-hand members the com- 
ponents of Dj(U), where U has for its components — M.vry, Mravr, 0. 

The last three have for their right-hand members the components of D^H), where 
H has for its components — vr%(mxz), —^(myz), M#W. And the six equations are at 
once obtained by applying the elementary formulae of section 18. 

60. It is, however, in the solution of problems, far better to avoid using those six 
equations, and simply to remember that the body's momentum U is Mwr in the direction 
of the velocity of the centre of gravity (f being its distance from the axis), and that the 
body's momentum-couple H has for its components — ro-2(»m), —sr%(myz), and M#W. 
Then the complete differential coefficients of U and H can be found at once according to 
the ordinary rules ; and those complete differential coefficients are by D'Alembert's 
principle respectively identical with the force and the axis of the couple to which the 
forces acting on the body may be reduced. 

Take for example the following well-known problem : — 

" Under what circumstances will there be no pressure on the fixed axis, supposing no 
external forces to act on the body % " 

Since there are no external forces nor pressures which act on the body, it follows that 
D,(U) and D,(H) must each equal zero. Therefore U and H are lines of constant mag- 
nitude and direction. Now the direction of U is that of the velocity of the centre of 
gravity, and would therefore vary, unless the centre of gravity were at rest. Hence the 
first condition is that the fixed axis passes through the centre of gravity. 

Again, since H is a line of constant length and direction, its components along and 
perpendicular to the fixed axis Oz must be lines of constant length and direction. 
Hence h z =M.k' i vT must be constant. Therefore ■m is constant, or the body revolves with 
uniform angular velocity. 

Moreover the components of H perpendicular to Oz are h x = —■m'%(mxz),h, J = —rwX(myz); 
and we have just seen that the resultant of these two components must be a line of 
constant length and direction. But as vs is constant, it is clear that that resultant has 
always the same components along the variable axes of x and y, and would therefore 
move with the latter, unless those components were always zero. Hence the second 
condition is that S(m#z)=0, %(myz)=0 ; in other words, the fixed axis must be a prin- 
cipal axis. It is evident also that the two conditions are sufficient, for they make H and 
U constant lines, and therefore they make D t (H) and D,(U) vanish, and consequently, 
by D'Alembert's principle, there are no forces acting on the body. 
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61. We now come to the case of a body moving about a fixed point O. Let m be, at 
time t, the angular velocity about the instantaneous axis, let A, B, C be the moments of 
inertia about the principal axes at O, and let vx x , w y , w z be the components of ^ along 
these axes. 

H, the axis of the body's momentum-couple about O, has, we have already seen, for 

its components 

h x =Avr x , h y =~E>Tff y , h^O&r 

If, then, we denote by Q the instantaneous axis, we know, from the fundamental propo- 

dJlx II tn Or ^s 
If l( to u *' Tt 



sitions in sections 35 and 36, that D,(H) is equivalent to ~^ || to Ox, -f || to Oy, 



~ || to Oz, together with the determinant of Q to H, Q denoting the instantaneous axis ; 

and moreover, that this determinant has for its components 

h z m v —\w z , or (C— B)^^. parallel to Ox; 
h x -ar z —h x Tn x , or {A.—C)vs x m z parallel to Oy ; 
h y nr x —h x 7!J r or (B— A)sya\„ parallel to Oz. 

Therefore the components of D ( (H) are 

C^+(B-A) V , 

But by D'Alembeet's principle D,(H) is the same as G, the axis of the couple result- 
ing from the external forces. If, therefore, L, M, N be the components of G, L, M, N 
must be respectively equal to the components of D,(H). Hence we have 

L=A^+(C-B)^, 

N=C^+(B-A>,^ 

We thus see that these well-known equations of Eulee are found at once by resolving 
D ( (H) along the principal axes, where H is the axis of the body's momentum-couple 
and has Ac7„ B^, Ota z for its components, and that they merely express the fact that G, 
the resultant of L, M, N, is identical with D,(H). 

62. The theory of the motion of a body about a fixed point can be more simply inves- 
tigated, and the problems connected with that theory can generally be more easily solved, 
by merely bearing in mind that G=D,(H) than by using Euleb's equations, which 
merely express that fact in one 'particular form ; for that form is not always the most 
convenient form, and is in all cases apt to conceal the fact which it embodies. 
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Take, for instance, the problem of a body rotating about a fixed point, no external 
forces acting on it. Here D,(H)=0, therefore H is a line of constant length and direc- 
tion. The motion of the body must therefore entirely depend upon the fact that, whilst 
the body moves about the principal axes with angular velocities w x , w y , w z , the line H, whose 
projections on those axes are respectively A^, IW y , Czzr z , remains throughout the body's 
motion the same in magnitude and direction. 

The length of H is evidently ^A VJ+BV^+C 2 ^. But the length of H is constant, 
say equal to h. Therefore 

AX+Bv;+ey:=tf. ........ (i.) 

Moreover, from section 35, we see that D«(H) is equivalent to A ~ \\ to Ox, 

B ^p || to Oy, C ^ || to Oz, together with det (Q, H) ; and since the last line det (Q, H) 

is perpendicular on the instantaneous axis 12, it follows that the resolved part of D^(H) 
along the instantaneous axis equals 

■st dt "S3- dt •& dt ' 

But this must equal zero, since D ( (H) equals zero, and since consequently its resolved 
part along any line is zero. Hence we have 

Therefore 

Ara^+Bwy+Cw* is constant, equal, say, to ¥. .... (II.) 

We have already seen that H is a line of constant direction ; and since its direction- 
cosines are proportional to A^, Bsr,, Cm e , it follows that the plane whose moment has 
direction-cosines which are proportional to the last three quantities is a fixed plane. 
This plane is the invariable plane. From this fact and the two equations (I.) and (II.), 
Poinsot's celebrated illustration of the motion of a body which rotates about a fixed 
point may be easily deduced in the ordinary manner ; but it is unnecessary to discuss 
the problem further, as it must be already sufficiently apparent that the body's motion 
entirely depends upon the fact that H is a line of fixed length and direction. 

63. On looking at Eulek's equations, we find that, when A=B=C, they take the 
simple form 

T A dlffx 

M=B ip, 

dt ' 

dt 
Moreover, when only A=B, then the third of Euler's equations becomes N=0~?. 

dt 

It may be interesting to trace the real meaning of these results. 
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Let, as before, G be the resultant of L, M, N, and let H be the axis of the body's 
momentum-couple about the fixed point, and let Q represent the instantaneous axis. 

Then G=L\(H), and D,(H) is equivalent to A ^ || to Ox, B^ || to Oy, C *£ [| to Oz, 

(It HZ Clt 

together with det (Q, H). 

Now, if A=B=C, then H, which has for its components A& x , Bsr y , CW,, evidently coin- 
cides in direction with Q, which has for its components tsr,, nr p ra> Therefore it follows, 
from the very definition of a determinant, that det (Q, H)=0. It is therefore because 

det (Q, H)=0 when A=B=C, that the components of D/H) are simply A *£, B ^, 

ut (Xt 

C ^~, and that Euler's equations take so simple a form. 

Gib 

Secondly, suppose only A=B. It is evident, from what has been just said, that 
N=C -^?+the resolved part along Oz of det (Q, H). Now the equation to the line H is 

Git 

x y z 



Awj Bit, Gsr, 

If then A=B, the projection of H on the plane of xy evidently coincides with the 
projection of Q on that plane. Therefore the lines H, Q, and the axis of z lie in the 
same plane. Hence it follows that the line det (Q, H), which by definition is perpendi- 
cular on Q and on H, is also perpendicular on the axis of z, and has therefore no com- 
ponent along that axis. We thus see that the reason why N=C ~, when A=B, is 
that in that case det (Q,, H) is perpendicular to the axis of z. 

64. In those cases in which there are more sets than one of principal axes at the fixed 
point, it is sometimes convenient to take moments about a set of principal axes, which 
are not fixed in the body. 

There is no difficulty in applying the same method to such cases. Let vr x , zs y , vt z be 
the angular velocities of the body about the principal axes Ox, Oy, Oz, and suppose those 
axes not to move with the body as if rigidly connected with it, but to move at time t about 
an instantaneous axis Q' with an angular velocity equal to the length of Q', and let the 
components of QI along the principal axes be w' x , vr^ •&',.. 

Let H, as before, represent the body's momentum-couple. We have seen that H has 

for its components 

h x =Am x , h^Biffy, h z =Cvt z . 

Therefore, according to the fundamental proposition in section 35, D/H) is equivalent 

to j. (Ato-.,) I] to Ox, j { (Bcr,,) || to Oy, j t ((V e ) || to Oz, together with det(Q', H); and looking 

at the formulae of section 36, we see that det(Q', H) has for its components 

Cos-^— Bw^ parallel to Ox, 
AsrX - CWgV, parallel to Oy, 
Bro-j,^— Aro^ parallel to Oz. 
MDCCCLXII. 3 Y 
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Therefore, if L be the moments of inertia of the forces about the principal axes, we 
have by D'Alembekt's principle, 

It is clear that, when there is more than one set of principal axes at the fixed point, 
either all three or at least two of the quantities A, B, C must be equal to one another. 
Suppose then A=B, then the axis of z, Oz is fixed in the body, and therefore m' x and w^ 
are clearly the same as ^ and ro- r And the last equations, therefore, become 

L=A^+^(<X-A<), 
M=A^+^(A<~CV 2 ), 

dt 
Ifweput^rr^+^y, the above equations become the same as those which are given 

in Routh's ' Dynamics,' page 134, where they are deduced from Eulee's equations. 

65. To the above equations, however, the same remark applies as has already been 
made with regard to Euler's equations. They merely express the fact that D 4 (H) has 
L, M, N for its components ; and it is far better in most problems to start with that 
simple fact, and, without using those equations, to choose any axes which the nature of 
the problem may suggest. 

Take, for instance, the problem of the top spinning upon a per- 
fectly rough plane. 

Let O be the fixed point, g the top's centre of gravity. Take O g 
as axis of z. Draw O a vertically, and take as axis of x a line per- 
pendicular to O z and in the plane z O a, and take as axis of y a line 
perpendicular on the plane z Ox, and therefore perpendicular on O g. 
The axes of coordinates are evidently principal axes. 

The components of H, which determines the body's momentum-couple, are A*r x , Bvr p 

Cvf z , vt x , vs y , w z being the angular velocities about the principal axes, and A being the 

moment of inertia about x and about O y, and C being the moment of inertia about O g. 

We have chosen O g so as to be perpendicular on plane a O x ; consequently the resolved 

part of H along O a is the sum of the resolved parts of A& x , CV,., and equals therefore, if 

we denote angle a O z by 6, 

— Asr x sin 6-\- Qrs z cos 6. 

Moreover, since O a has a fixed direction, the differential coefficient of the last expres- 
sion is evidently the resolved part along O a of D,(H). But this must equal zero by 
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D'Alembert's principle, since there are no forces acting on the top which have any 
moments about the vertical. Hence we have 

d 

j t (Czif g cosd— ■Aa.sia.fy—O (I.) 

Therefore 

Oa t cos 6—Azs x sin 6 is constant. 

Moreover, since two of the principal moments of inertia are equal to one another, it 
follows from section 58 that the sum of the moments of the external forces about Oz, 

the axis of unequal moment of inertia, is equal to C -^?. Hence, as the forces have no 

(hi 

J 

moment about 0#, we have C —f=0. Therefore & e is constant, This relation, together 

with the equation (I.) and the equation of vis viva, solve the problem. We, namely, 
obtain the three equations 

Oa x cos 8— Avs x sin &=h, 

Aw£+ Affy+Cnr^= —2gi cos 0+c, 
where b=Og, and c is some constant. 

66. In some few cases it may be convenient to take moments about lines which are 
fixed in the body but which are not principal axes. Let then H have for its components 
along the rectangular axes of coordinates h x , h y , h s respectively, then we know, from 
section 53, that 

h x = Aa x — B'ray— C'jsr y , where A' = %(mys), 

h y =Bzs y — Ozfs—A'zfg, B'=2(»m), 

h z =.Ofi e — A'w v — B'tf,,, G=X(myx). 

Therefore, if L, M, N be the moments of the forces about the axes, we have, as before, 

If, on the other hand, the axes of Ox, Oy, Oz are not fixed in the body, but rotate with 
an angular velocity whose components are <a m is v , &'„ then we have, in a similar manner, 

L= ^*+hj*,- V- &c - 

The above equations are somewhat more general than those given by Liouville in 
his Journal of 1858, and are, as we have just seen, at once obtained by applying the 
fundamental formulae of section 36. 

67. I will now show how the principle of vis viva may be easily proved for a body 

3t2 
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moving about a fixed point without assuming the principle of virtual velocities, and is in 
fact a very simple deduction from D'Alembeet's principle. ^ B 



■4 



» 



Let O i be the instantaneous axis, about which the body is rotating 
at time t with an angular velocity vr. Let m be a particle of the body, 
let its mass be to, its distance from O i, r, and its velocity v in the 
direction to n. 

Suppose now a force P to act at to. Then, since to n is perpendicular on the plane 
i O to, it is clear from statics that the moment of P about O * is the moment of the 
resolved part of P along to n. It is therefore, if <p be the angle which the direction of P 

V 

makes with to n, equal to rP cos <p, or - P cos <p. 

Suppose then P to be the moving force of the particle to. Then its resolved part 

dv ~ dv 

along the velocity to n is of course m-j t , so that P cos <p=m-j t , and therefore the moment 

v dv 

of the moving force about i equals to - -r.. Consequently the sum of the moments of 

the moving forces about O i equals - % (mv-jA. But this sum, by D'Alembeet's prin- 
ciple, equals the sum of the moments of the external forces about (H. Now we have 
already seen that the moment of any force P about i is - P cos <p, so that, if P represent 
an external force acting on the body, the sum of the moments about O i of the forces 
acting on the body equals - 2 (Pa cos <p). Hence we have 

Therefore 

S(m<) =2§dt%(Pvcos<p). 

This equation embodies the principle of vis viva ; for it is evident that, if the compo- 
nents of P be X, Y, Z, then 

P*cosp=X^+Yj+Z 7 ,. 

68. The same result may also be obtained by analysis; and it may be worth while to 
notice that each step in the analytical proof is exactly equivalent to the corresponding 
step in the above geometrical proof. This correspondence between the steps in analytical 
and geometrical demonstrations is one of the most striking features of modern analytical 
geometry, and would, as we have already attempted to show, present itself generally in 
analytical mechanics, if more attention were paid to the interpretation of the equations 
and formulae which are employed. 

The analytical proof is as follows : — 

Let X, Y, Z be the components of any one of the forces acting on the body, and sup- 
pose that force to act at a point (a?, y, z) of mass to. Let or,, &„, w s be the angular velocities 
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of rotation about the axes of coordinates, which are here supposed to be fixed in space. 
Then it is clear that 

dx dy dz 

■^=zs$Ts y —yv, z , jj—XT3 z —Svr x , j t =ynr x —X7ff r 

Therefore 

may be put into the form 

2(Xs-Z*K+:S(Ytf-X#H+S(Zy-Y;^ (I.) 

But, by D'Alembbet's principle, 

S(Xz— Zx)—%m (jgj z—^x), 
2(Yx-Xy)=2m(^ X -§y), 

S(Zy-T*)=a»(^-g[.). 

Substituting these expressions in (I.), we obtain 

_ /«P* d?z \ , ^ /d 2 y <fe 2 \ , - /(Pz d*y \ 

^^ m {d¥ z -d¥ x )+ u ^ m {d¥ x -d¥y)+^^ m \dt^-d¥ z )^ 

which again can be put into the form 

(d 2 & d?y d?z 1 

S»»|^-(«w f — yw^+^Ka*,— «w.)+^s(yw.— aw,)j, 

/ <Px dx d^y dy d*z dz\ 
= * m \dl?d7 + d¥d7'T~~dtJt)- 

Hence it follows that 



and therefore 



2 ^ X ¥+ Y <ft+ Z dt)~ 2 ' m \df dt^dt* dt+dPdt)* 
re 

M(ty+ (ty+ (lyHj^+^+^j' 



which is the equation of m's wya. 

It may be observed that the same proof may be quite easily extended to a body moving 
freely, by decomposing the original motion into a motion of translation with the velocity 
of the centre of gravity and a rotation about the centre of gravity. 

69. The two proofs which have just been given of the principle of vis viva are both 
founded on the fact that the sum of the moments of the moving forces of the body's 

particles about the instantaneous axis is equal to - 2 ( mv jX or to ^ j ( 2 (mi; 2 ). This 
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fact is the reason why the equations of vis viva can he obtained by multiplying Euler's 
three equations byro^, vr y , m z respectively, and by adding together the products so 
obtained ; for in performing those operations we are in reality finding the sum of the 
moments of the forces about the instantaneous axis whose direction-cosines are 

^f S fl* 

•or ' ts ' w * 

70. If Gi be the sum of the moments of the impressed forces about the instantaneous 
axis, then by D'Alembbet's principle G x equals the sum of the moments of the moving 
forces about the instantaneous axis. Hence it follows from the preceding section that 

This equation is often useful. 

For instance, since 2(mfl 2 )=2(msrV), where r is the distance of a particle from the 
instantaneous axis, it follows that 2(ww> 2 )=Isr 2 if I denote the body's moment of inertia 
about the instantaneous axis. Therefore 

_<fe vrdl 
— ~dt~T'2Tf 

Now, if the instantaneous axis were fixed in space, we should evidently have 

Therefore the only cases in which we can take moments about the instantaneous axis 

as if it were fixed in space are when ^ "^=0, or when the moment of inertia about the 

instantaneous axis is constant. This proposition is useful in solving problems concerning 
rolling cones, and is usually deduced by analysis from Euler's equations. 

71. I will give two more examples of the advantage of the method I have employed 
in these pages. 

Let & x , Wy, vr z be the angular velocities of rotation of a body about three rectangular 
axes which are fixed in and move with the body, and let a, b, c be the direction-cosines 
with respect to those axes of a line which is fixed in space. Take on the latter line a 
point P at a unit of distance from the origin. The velocity of the fixed point P is zero. 
Now, as its components along the moving axes are a, h, c respectively, it follows from one 
of our elementary propositions that 

da , 

equals the component of P's velocity along Oir, and therefore equals zero. 



AJSTD DETEEMINANTS OP LINES. 



507 



Hence 



Jt=^- Cm p 



a formula which is generally deduced as the result of somewhat long analytical work. 

72. Secondly, in order to give a striking example of the manner in which the theory 
of the determinants of lines explains and shortens analytical processes, I will give the 
following direct proof of Eulek's equations. 

Take the principal axes at the fixed point as the axes of coordinates. Let %, v y , v z be 
the components of the velocity V, andy^, f p f z those of the acceleration F of a particle 
(a?, y, z) of mass m. 

We have seen that the fact of the acceleration being the complete differential coefii- 
cient of the velocity leads at once to the three following equations : — 

j, dv x , 

f*=it +v * 1ff ~ v °™ ! °> \ • * • • P') 

n dVs, . 

Putting then for brevity's sake f' x for v z vTy—Vy-m z ,fy for v x ™ z —v z m x ,f' z for v y ™ x —v x rss y , 
we have 



(2.) 



n __dv x , _pi \ 

> x ~ dt" rj x ' 

Jy— dt -TJ r 

f dv z ,j;i 

Jz —~dt' rji - 

Now the sum of the moments of the moving forces about Ox equals 

2m(/^-^)=2m(^2/-J^)+2m(/^-/» .... (3.) 

Let us first investigate the expression %m(~p y—^ #)• Evidently 

v z =yvr x —xis s , and » } =aw e -2w, 

And as the axes m'@ve with the body, w, y, z do not vary with the time, and we there- 
fore obtain 



Therefore 



fa fes__ fey dv v dts z _ d-st x 

dt ~~ y dt dt' dt dt dt' 
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But as the axes of coordinate axes are principal axes, %(mxy)=0, %(mxz)=Q. 
Therefore 

A being the moment of inertia about Ox. 

We have still to investigate the expression %m(f^.y—f 9 z) where ^'=^5^— v t m s , and 

We obtain by substitution 

%m(fly—f y z)=mXH v yy+ v ^)—^M mv ^)—^^{ mv ^)- 
But since 

v„=z^ 9 — yet,, 

^=W 2 — Zvr x , 

it is evident that for principal axes we have 

%(mv s y)=0, X(mv z z)=z0, 

t(mv^)= —a,g(mf), X(mv x z)=rs^(mz s ). 

=w 3 ,Trr ;! .(2m(y+* 2 )— Sw(0 2 4-* 2 )) 
=^^(C-B), 

where C and B are the moments of inertia about the axes of z and y. 

Substituting then this last expression and the expression in (4.) in equation (3.), we see 
that the sum of the moments of the moving forces about Ox equals 

Hence by D'Albmbekt's principle, if L, M, N be the sum of the moments of the 
impressed forces about the axes of x, y, 2, we have 

L=A^+(C-BK^. 
Similarly, 

m=bJm-(a-ck^, 

N=C §5-HB-AK*v 

I will now show the full import of each of the steps in the above analytical proof. 

We have seen, in section 39 of Chapter III., that the acceleration of the particle is 
the result of the acceleration det (D t (Q), K), and the acceleration det (Q, V). If then 
the particle's mass be m, its moving force is represented by 

m det (D,(Q), E)+«i det (Q, V). 
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Now the above analytical proof merely shows that the sum of the moments about 
the coordinate axes, of the moving forces of which m det (D t £2, R) is the type, are 

respectively A — , B -^, C ~, and that the sum of the moments about the coordinate 
r J dt' dt' dt' 

axes of the moving forces of which m det (Q, V) is the type are respectively (C — Bjsr^, 

(A— C>*a, (B— A)*,*,. 

73. I will next proceed to show how these results can be obtained far more briefly by 

applying the propositions concerning the determinants of lines. In the first place, if we 

put, for D t (Q), P, and suppose P to have for its components^,, p $ , p z parallel to the axes, 

then m det (D^Q), R) has for its components parallel to z and to y, 

mfyp x —xp t ) and m(xp z —zp x ), 
and therefore the moment of m det (D t (Q,), R) about the axis of x equals 

m \ (yP°> — X P^—( X P« — S P^) Z \ =m{z i J rf)Px—my%P ll — mxzp,,. 
Therefore, if we take the sum of these for all the particles and remember that the 
axes are principal axes, that sum will equal p x 2m(z 2 -\-y 2 ). Now we have already 

proved that p x , the component parallel to the axis of x of D,(£2),= -£. Hence we see 

that the sum of the moments about Ox of the moving forces of which m det (D t (Q,), R) 

is the type is A -j~, and that this follows from the properties of principal axes, and from 

the fact that the component of D,,(Q) parallel to the axis of x is -£. 

In the second place, we have already seen in section 30 of Chapter II. that the 
moment about the axis of x of det (Q, V) equals — rav x cos <p, v x being the component of 
V, and <p the angle between the radius vector and V. 

But —rixv x cos <p equals evidently — ^(tf^+^+sw,), and v x equals ZTX y —yw z . If 
then we observe that '2(mxv x )=0, H,(myv x )=—'2(my i )w ! ., '%(mzv x )='2(mz 2 )& !/ , it follows 
easily from the above that the sum of the moments about Ox of the moving forces, of 
which m det (Q, V) is the type, equals 2(m^ 3 ) sysr z — 1,(my 2 ) w^ y :=(C— B)^^. This 
last proposition may be also proved in a different manner, which will show its connexion 
with the proof first given of Euler's equations. 

The moment-axis, with respect to the origin, of the acceleration det (Q, V) is 
det {R, det (Q, V)}, which, as we see from section 30 of Chapter II., equals a line oppo- 
site to V and of length raw cos <p ; but it follows from the same section that, since Q is 
perpendicular on V, det {Q, det (R, V)} equals a line opposite to V and of length 
wvr cos <p. Hence we have 

det {R, det (Q, V)}=det {Q, det (R, V)}. 

Let then % denote the operation of taking the complete sum of lines. Then it follows 
from the last equation that 

2m det{R, det (Q, V)} =%m det{Q, det (R, V)} 
= det(Q,2mdet(R,V)). 

MDCCCLXII. 3 Z 
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Now 2m det (E, V) equals the complete sum of the moment-axes of the momenta, or 
is, in other words, the axis of the body's momentum couple H whose components are 
AjTj,, Bsr y , Cw 2 . Hence we see that 

Xm det (R, det (Q, V))=det (Q, H), 
and therefore the sum of the moments about Ox of the moving force, of which 
m det (Q, V) is the type, is the component parallel to x of det (12, H), and equals 
therefore Cnr^— Bay^=(C— B)^w 2 . 



